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We extend chiral perturbation theory to study a meson gas 
out of thermal equilibrium. We assume that the system is 
initially in thermal equilibrium at a temperature Ti < T c and 
work within the Schwinger-Keldysh contour technique. To 
lowest order and in the chiral limit, we use a nonlinear sigma 
model with a time-dependent pion decay function f-n(t), which 
we consistently define in terms of the axial current two-point 
correlator. As a first application, we analyse up to one 

loop order and find the T% dependence of the lowest coefficients 
in its short-time expansion. We discuss the applicability of 
our results to the time evolution of the plasma formed after a 
heavy-ion collision and discuss the interpretation of our model 
in a curved space-time background. 



I. INTRODUCTION 

The chiral phase transition plays a fundamental role in 
the description of the plasma formed after a relativistic 
heavy- ion collision (RHIC). For Nf massless quarks (the 
chiral limit), QCD is invariant under the chiral group 
SUL(Nf) x SUft(Nf). That symmetry is spontaneously 
broken at zero temperature to the vector group SUv{Nt) 
(isospin) , the Nambu-Goldstone bosons (NGB) being the 
lightest mesons (pions for Nf — 2 plus kaons and eta for 
Nj = 3), which are the relevant low-energy degrees of 
freedom. The light quark masses are introduced pertur- 
batively. This picture has proven to be very successful 
in describing hadron observables (like decay constants or 
scattering amplitudes) , for low energies and small masses, 
through the Goldstone theorem and current algebra. On 
the other hand, QCD at this scale is strongly coupled 
and, therefore, if we want to describe properly the meson 
dynamics we need an effective low-energy model. Since 
SU(2) x SU{2) -► SU{2) is isomorphic to 0(4) -► 0(3) 
then for Nf = 2 a candidate for the effective theory is the 
0(4) linear sigma model (LSM), whose classical poten- 
tial has the familiar Mexican hat form. Its fundamental 
fields are the three pions and the a, which acquires a 
nonzero vacuum expectation value (VEV) v. However, 
one drawback of the LSM is that its coupling constant 
A becomes large at low energies, so that perturbation 



theory in A is meaningless and alternative perturbative 
expansions must be considered, such as taking a large 
number TV of NGB. Besides, the LSM does no longer have 
the QCD chiral symmetry breaking pattern for Nf = 3. 
An alternative approach is to construct an effective the- 
ory as an infinite sum of terms with increasing number 
of derivatives and with all the QCD symmetries, only for 
the NGB, which are included in a SU(Nf)-valued field, 
since SU(Nf) is the coset space of the QCD chiral sym- 
metry breaking pattern. To lowest order, such theory 
is the so called nonlinear sigma model (NLSM), which 
allows, not only to rederive the current algebra predic- 
tions straightforwardly, but to obtain corrections to them 
in a consistent framework. Such corrections come both 
from NGB loops and higher order lagrangians and can be 
renormalised order by order in energies, yielding unam- 
biguous finite predictions for the observables. This con- 
stitutes the so called chiral perturbation theory (ChPT) 
(llpl where the perturbative expansion is carried out in 
terms of the ratio of the meson energy scales of the the- 
ory (masses, external momenta, temperature and so on) 
and the chiral scale A x ~ 1 GeV. For a review of low- 
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energy QCD and effective lagrangians we refer to 
and references therein. 

In thermal equilibrium at finite temperature T, the 
chiral symmetry is believed to be restored at T c ~ 150- 
200 MeV in a phase transition, which is very likely to 
be of second order for Nf = 2 It should be stressed 
that ChPT cannot strictly establish the existence of a 
phase transition since it provides an expansion in T 2 /A x 
for the order parameter VT. Nonetheless, it predicts the 
correct behaviour of the observables as T approaches T c 
(see below) and allows to study the meson gas thermo- 
dynamics at low temperatures [||J^] . To describe the sys- 
tem near T c , the LSM seems to be a better choice, since 
it undergoes a second-order phase transition in the mean 
field approximation, v(T) becoming the order parameter. 
However, T c in the mean field LSM is independent of A, 
which suggests that the NLSM should be equally valid 
to study the phase transition, which is indeed the case 
provided one works in the large N limit S . 

In the chiral limit and in equilibrium at temperature 
T, the only parameter of the NLSM lagrangian is a T- 
independent constant / with dimensions of energy, re- 
lated to the physical temperature-dependent pion decay 
constant as U(T) = f(l + <D(p 2 )), where f w ~ 93 MeV 
at T = (A x ~ 471"/^, which is the typical loop factor) is 
customarily measured in leptonic decays ir — > lv[. Here, p 
denotes generically any pion momentum scale, including 
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temperature and by 0(p) we really mean 0(p/A x ), as 
we explained before. Besides, every pion loop counts as 
0(p 2 ) § and all the infinities coming from them can be 
absorbed in the coefficients of higher order lagrangians, 
which is also true for T ^ . The next to leading or- 
der (NLO) one-loop corrections induce the temperature 
dependence of f_rr{T) and for Nf = 2 in the chiral limit 
the result is H[| 



fl{T) = f 1 



6p 



(1) 



A few remarks are in order here. Firstly, at T ^ 0, due 
to the loss of Lorentz covariance in the thermal bath, one 
should actually consider two independent pion decay con- 
stants, /* (temporal) and (spatial), which in general 
will be complex, their real and imaginary parts being 
related respectively with the pion velocity and damping 
rate in the thermal bath 0). However, to NLO, one has 
/*=/* = U(T) in (@). Notice that (@) suggests a criti- 
cal behaviour at T c ~ V6fn, since f^{T) behaves roughly 
like v(T), although they are different objects Despite 
(|I|) being just the lowest order in the low temperature ex- 
pansion (and thus not being too accurate to extrapolate 
it up to T c ) it predicts the right behaviour and a rea- 
sonable estimate for T c . Another important comment is 
that one has to be careful when defining f n in a medium, 
trying to extend naively low-energy theorems to T ^ 
(see section [V). Nevertheless, there is a consistent way 
of defining f n in the thermal bath from the axial-axial 
correlator spectral function |^] as it will be explained in 
section 



IV 



The equilibrium assumption is not realistic if one is in- 
terested in the dynamics of the expanding plasma formed 
after a RHIC, where there are several observable non- 
equilibrium effects. One of them is the formation of the 
so called disoriented chiral condensates (DCC), regions 
in which the chiral field is correlated and has nonzero 
components in the pion direction |Tc[ |. As the plasma 
expands, long- wavelength pion modes can develop insta- 
bilities that grow very fast as the field relaxes to the 
ground state. If that picture is correct, there would be 
observable effects, such as coherent pion emission [ p"T| |. 
The question is whether those domains can grow fast 
and large enough to emit a significant number of pions. 
A necessary, but not sufficient, condition is that the sys- 
tem is out of thermal equilibrium. This problem has 
been extensively studied in the literature, using various 
non-perturbative and numerical approaches within the 
LSM p3|— |Tq] . In all these works, thermal equilibrium at 
Ti > T c is assumed for t < and then some mechanism 
drives the theory off equilibrium for t > 0, encoded in 
the time dependence of the different lagrangian parame- 
ters. The simplest choice is an external quench changing 
the sign of the mass squared m 2 of the 0(4) field so that 
modes with k 2 < m? become unstable pl[-Jl4]|. A more 
realistic approach is to describe the plasma expansion in 
proper time and rapidity coordinates, which is analogous 



to study the system in a certain space-time background 
metric, assuming either boost [ |l5| or scale 16 invariant 
kinematics (cylindrical and spherical expansions, respec- 
tively). During the time evolution, the pion modes prop- 
agate as if they had a negative mass squared (unstable 
long- wavelength modes), therefore yielding exponential 
growth with time of the pion propagator. The general 
conclusion is that the typical size of the formed DCC 
regions is at most 3 fm, which appears to be not big 
enough to contain a large observable number of pions, 
and the time during which the relevant plasma expan- 
sion takes place varies between 5-10 fm/c ]15]Jl6]] . An- 
other important nonequilibrium effect is the photon and 
dilepton production Jl7],[l8| . As these particles are weakly 
interacting, they decouple from the hot plasma carrying 
important information about the dynamics, so that their 
final spectra, dominated by meson decays, is a promising 
signal to analyze nonequilibrium aspects [ p~7| . A recent 
suggestion is that a significant number of those photons 
might have been produced by 7r° decay during a typical 
DCC nonequilibrium evolution Jl^ ]. 

Our basic objective in this work will be to construct 
an effective ChPT-based model to describe a meson gas 
out of thermal equilibrium, as an alternative to the LSM 
approach. Our only degrees of freedom will be then the 
NGB and we will consider the most general low-energy 
lagrangian compatible with the QCD symmetries. We 
will restrict here to Nf — 2 and to the chiral limit. This 
is the simplest approximation we can consider and it al- 
lows to build the model in terms of exact chiral symmetry. 
The lowest order lagrangian is then the NLSM. In princi- 
ple, we will be interested in the physical regime in which 
the system is not far from equilibrium (we will be more 
specific about this point below), where an expansion in 
derivatives is consistent. One of the novelties of our ap- 
proach is to exploit the analogy between near-equilibrium 
systems and ChPT. 

The structure of this work is the following: in section 
H we introduce our non-equilibrium NLSM and estab- 
lish the non-equilibrium ChPT, discussing how the dif- 
ferent symmetries are realised. We also discuss the in- 
terpretation of the model as a curved space-time QFT, 
paying special attention to rcnormalisation. The lead- 
ing order non-equilibrium pion propagator and Lehman 
spectral function are also analysed in that section, which 
also contains most of our notation and conventions. In 
section III, we derive the one- loop NLO correction to the 
propagator, which will be needed later. In section |^ 
we provide a consistent definition of the nonequilibrium 
time-dependent pion decay functions (PDF) extending 
the equilibrium pion decay constants and we derive them 
up to one loop in ChPT. One of the motivations to con- 
centrate on /?!■(£) first is that to one loop in ChPT we 
expect that all physical observables are obtained from 
the tree level ones by replacing f(t) — > fn(t), as it hap- 
pens indeed in equilibrium. We analyse the short-time 
behaviour of f n (t) and estimate the relevant time scales 
involved, discussing the role of unstable modes. The issue 
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of axial current conservation is also analysed in section 
IV . Finally, in section 0, we present our conclusions and 
discuss some open questions, as well as further applica- 
tions of our approach. 



II. THE NLSM AND CHPT OUT OF 
EQUILIBRIUM 
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FIG. 1. The contour C in complex time t. The lines Ci 
and C2 run between ti + ie and tf + ie and tf — ie and ti — ie 
respectively, with e —* + . 

We will assume that the system is in thermal equilib- 
rium for t < at a temperature Tj < T c and for t > 
we let the lagrangian parameters be time-dependent (we 
are also assuming that the system is homogeneous and 
isotropic). The generating functional of the theory can 
then be formulated in the path integral formalism, by let- 
ting the time integrals run over the Schwinger-Keldysh 
contour C showed in Fig.[j] []l9|-^2||. We will eventually 
let ti — > — 00 and tf — > +00, although we will show be- 
low that our results are independent of U and tf. We 
remark that, even in that limit, the imaginary-time leg 
of the contour has to be kept, since it encodes the KMS 
equilibrium boundary conditions ]20| , p^ |. 

With the above assumptions, our low-energy model 
will be the NLSM with / becoming a real function f(t), 

S[U] = [ d 4 x tM. tr d„tf(x,t)d»U(x,t) (2) 
Jc 4 

with J c d 4 x ee J c dtJ d 3 x, f(t < 0) = / and U(x,t) € 
SU{2) is the NGB field, satisfying U(x, ti+ifc) = U(x, U) 
with Pi = T f . We will parametrise U as 

U(x, t) = -ffi{ [/ 2 W - ^ V2 1 + ^ a } (3) 

where ir 2 = 7r a 7r a , ir a (x,t) the pion fields satisfying 
ir a (ti + ifli) = ir a (ti) (a — 1,2,3) and / and r a are the 
identity and Pauli matrices. In terms of the 7r° fields, the 
action (0) becomes 



S[n] = J d^x X - (d^ a d^TT 



P(t)-^ 



/*(*)- 7T* 



m 
at) 



—1 rx a-o 
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(4) 



Note that by including f(t) in the parametrisation of 
the NGB field we recover the canonical kinetic term in 
the action. Other choices amount to a time-dependent 
normalisation of the pion fields, which should not have 
any effect on the physical observables (see section [Tv| ). 
For instance, we could redefine 7r a = ir a f (0) / f (t) , so 
that for the Tr a (x,t) fields the action is the equilibrium 
NLSM multiplied by the time-dependent scale factor 
f 2 (t)/f 2 (0), which is directly interpreted in terms of a 
curved space-time background, as we will see below. 

The action (^|) is manifestly chiral invariant, i.e, under 
U — * LUW with L and R arbitrary constant SU(2) ma- 
trices. Notice that, as we are working in the chiral limit 
by assumption, we have not considered a pion mass term 
in the action. Such term would break explicitly the chi- 
ral symmetry, being still invariant under vector isospin 
L = R transformations for m u = nid- The conserved 
vector V!f and axial-vector Af. Noether currents for the 
chiral symmetry can be derived by gauging the symme- 
try and then taking functional derivatives with respect 
to the external sources . Extending that procedure to 
our non-equilibrium action (^J) we get 



A«(x, t) = i^-tx [r a (U^d^U - Ud^)] 
y;(.f, t) = -i^P-tr [r a (U^d^U + Udprf)] 



(5) 



The model is also parity invariant. Lorentz covariance 
is lost, because we have chosen the frame in which the 
thermal bath is at rest for t < 0. The time translation 
symmetry is also broken at noncquilibrium. 

Let us now discuss how to implement a consistent 
ChPT out of equilibrium. The new ingredient with re- 
spect to the equilibrium case, is the temporal variation 
of f(t). We will then consider 



m 
Pit) 



Oip), 



m [fit)? 



Pity p(t) 



Oip 2 ), 



(6) 



and so on, the rest of the chiral power counting being 
the same as in equilibrium, i.e, Tj = 0(p), dn = 0(p 2 ) 
and every pion loop is 0(p 2 ). Therefore, in our approach 
we treat the deviations of the system from equilibrium 
perturbatively, following the ChPT guidelines. Thus, 
with this chiral power counting, we can expand our ac- 
tion (H) to the relevant order in pion fields and calculate 
any observable by taking into account all the Feynman 
diagrams that contribute to that order. The loop di- 
vergences should be such that they can be absorbed in 
the coefficients of higher order lagrangians, which in gen- 
eral will require the introduction of new time-dependent 
counter-terms (see below). 

Notice that according to (^|), we can always describe 
the short-time non-equilibrium regime, just by Taylor ex- 
pand f{t) around t = 0. In fact, for times t < J" 1 , the 
Taylor expansion of f(t) is equivalent to a chiral expan- 
sion, since then f(0)t/f = 0(p), /(0)i 2 // = 0(p 2 ) and 
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so on. In case we are interested only in short times, the 
specific form of f(t) will not be important (see section 
IV), but we want to stress that the conditions (||) do not 
constrain us to work at short times, but just to remain 
close enough to equilibrium. 

To leading order in the pion fields, the action ([|) reads 

S [n] = - [ d 4 xl-Tr a (x, t) [□ + m 2 (t)] n a (x, t) (7) 
Jc 2 



[) fit) 



(8) 



where we have integrated by parts, using the boundary 
conditions for the pion fields. Thus, we see that the lead- 
ing order non-equilibrium effect of our model can be writ- 
ten as a time-dependent mass term for the pions. As 
we commented before, time-dependent field masses are 
a common (and welcomed) feature of out of equilibrium 
models [fl2|-fu| . In our approach, that mass parametrises 
the deviation of the plasma from equilibrium through 
the temporal variation of fit)- Notice that m 2 (t) can 
be negative, so that our model accommodates unstable 
pion modes, whose importance we have discussed before. 
The effect of those modes should not be important for 
the short-time evolution, but for longer times we expect 
them to influence significantly the dynamics. 

We remark that the axial current is classically con- 
served despite the existence of the time-dependent mass 
term. For instance, to leading order we have, from (^]), 



K = -fit)d^ a + W(*K + o(tt 3 ) 



(9) 



which satisfies d^A^ = using the equations of motion 
to the same order [□ + m 2 (t)] n a = 0. Had we included 
an explicit pion mass term m w ~ 140 MeV, the axial 
current would not have been conserved and the effective 
mass term in the leading order action would have been 
proportional to m 2 , + m 2 (t), so that the onset for insta- 
bilities would be rather m 2 (t) < — m^. 

We will now rephrase our model in the language of 
curved space-time background QFT pa ] , which will turn 
out to be a very useful analogy. For that purpose, let 
us consider the NLSM in an arbitrary background space- 
time metric g M „ p4| ], 

S g [X]= [ d^x\^-gg^ {d^n a d u n b (S ab 
Jc z 

+ /2( of^ 2 )}+^^] , (10) 

plus tt independent terms, where g is the metric deter- 
minant and the term £Sr[it, R] accounts for the possible 
couplings between the pion fields and the scalar curva- 
ture R(x) (like R(x)(j) 2 for a free scalar field <f> ||). To 
compare with our nonequilibrium model, we will choose 
a spatially flat Robertson- Walker (RW) space-time: 



with a(t) the scale factor. For t < we choose the 
space-time as Minkowskian, so that a(t < 0) = 1 and 
the system in thermal equilibrium at temperature T^. 
Then, if we change coordinates to the conformal time 
rj defined by t — a{rj')dri' , the metric becomes ds 2 — 
a 2 {r/)[drj 2 — dx 2 \ and we see that the action in (|Io| ) for the 
minimal coupling case i = is nothing but S[U] in (|), 
using the n parametrisation discussed before and iden- 
tifying a(t) = f(t)/f(0). That is, our non-equilibrium 
model in the chiral limit is equivalent to a RW spatially 
fiat space-time in conformal time, that starts expand- 
ing (or contracting) at t = with a rate given by the 
scale factor /(i)//(0) and couples minimally to the mat- 
ter fields. Our model is then not only suitable as an ef- 
fective theory for a RHIC environment, but also in a cos- 
mological framework. Nonequilibrium models for RHIC 
in which the plasma expansion is also parametrised in 
a background metric have been analysed in Jl5|,[l6| (see 
comments in section |). 

Let us now consider £ ^ 0. The lowest order Sr term 
we can construct has the form of an effective mass term, 



S R [7r,R]=R(t)f 2 (0) 



/ 2 (0) 



1/2 



(12) 



ds 2 = dt 2 



\t)da 



(11) 



so that it breaks explicitly the chiral symmetry. Let us 
then pick up the lowest order contribution in pion fields 
Sr = —Rn 2 /2 in (|l2|). Then recalling the expression 
of the scalar curvature for the spatially flat RW metric 
in conformal time R{ri) — 6/(?7)// 3 (?7), we see that the 
curvature term exactly cancels the m 2 (t) term in (Q) , af- 
ter rescaling n(x) — /(0)7r(cc)//(t), provided we choose 
£ = 1/6. But that value of £ is precisely the only one 
making the theory scale invariant |p3[ , so that in that 
case the leading order lagrangian would be scale invari- 
ant but not chiral invariant. Thus, the lagrangian chi- 
ral and conformal symmetries are not compatible in a 
curved background pH or, equivalently, at nonequilib- 
rium. Thus, for £ = we interpret the mass term m 2 (t) 
as the minimal coupling with the background, which for 
TO^r = yields chiral invariance of the lagrangian, though 
breaking the conformal symmetry. This is the choice we 
will adopt here, since we want to preserve chiral symme- 
try. 

A very important result is that all the one-loop diver- 
gences arising from (|l0|) can be absorbed in the coeffi- 
cients of the 0(p 4 ) lagrangian which includes new 
chiral-invariant couplings of pion fields with the curva- 
ture [p4f . In the chiral limit, for instance, there are 
two such new terms and therefore two more undeter- 
mined constants in addition to the Minkowski case ones 
(see |24) for details). Therefore, following this analogy, 
we know how to construct the higher order nonequilib- 
rium lagrangians that eventually will absorb all the di- 
vergences we may find in our analysis. We will come back 
to this point below. 

From the generating functional we can derive in the 
standard manner, the Green functions ordered in time 
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along the contour C 
the pion propagator 



22) . The two-point function defines 



G ab {x,y) = -i<T c n a (x)TT a (y) > 



(13) 



with Tc the time-ordered product along C. To leading 
order the pion propagator is G^ix, y) = S ab Go(x, y) (by 
isospin invariance), with Go(x,y) the solution of the dif- 
ferential equation 



{n x +m 2 (x°)}G (x,y)=-S^(x-y) 



(14) 



As for the boundary conditions, firstly we have to im- 
pose KMS equilibrium conditions on the imaginary leg 
of the contour, that is, 



Gq (x , U - i Pi ; y) = G$ (x, U ; y) 



(15) 



with the advanced and retarded propagators defined as 
customarily along the contour C . In addition, the prop- 
agator must be continuous and differentiable Vt € C, so 
that the solution is uniquely defined [f^oj . Therefore, in 
our case, we will demand that the solution of (|l4|) for 
t — > + matches that for t < 0, which is the equilibrium 
solution, given in momentum space by p5[ : 

Gg{g a ,w 9 )=(%{q ,w q ) [9 c (t - 1') + n B (q )] (16) 

where w 2 = \q\ 2 , ub{x) = [exp(/3x) — 1] _1 is the Bose- 
Einstein distribution function, 9q is the step function 
along the contour C and = — 27risgn(q )^('?o — w g)- 

We will make use of two different representations 
for the non-equilibrium Green functions. Notice that 
G(x, x') — G(t, t', x—x') because of the lack of time trans- 
lation invariance. Therefore, we will define, as customar- 
ily, the "fast" temporal variable t — t' and the "slow" one 
r ee (t + t')/2, so that F(q ,uj q ,T) and F(w q ,t,t') will 
denote, respectively, the fast and mixed (in which only 
the spatial coordinates are transformed) Fourier trans- 
forms of F(x,x'). The fast Fourier transform depends 
separately on go aud <-o q because of the loss of Lorentz co- 
variance and has the extra nonequilibrium r-dependence. 

In the mixed representation, ( fTi] ) becomes 



d 2 
df 2 



+ LU 2 + m 2 (t) 



G (cj q ,t,t') = -S{t-t') (17) 



The general solution of ( |l7j ) cannot be found analyti- 
cally for m 2 (t) arbitrary, but it can formally be written 
as a Schwinger-Dyson integral equation in terms of the 
equilibrium solution as 



G (u q ,t,t') = G e q (cj q ,t-t') 

+ / dzm 2 (z)Gl q {uj q ,t - z)Go(uj q ,z,t') 
Jc 



(18) 



with Gl q (uj qi t-t') the solution of (|17j) with m 2 {t) = 0. 

Another object of interest for our purposes is the 
Lehman spectral function, defined as pl|,p5t 



p(x,y) = G > (x,y) - G < (x,y) 



(19) 



Recall that in equilibrium and to leading order, the 
spectral function is p 1 ^ in (|l6|). Let us now dis- 
cuss some useful properties of the spectral function. 
Firstly, by definition G(x, y) — G{y, x) and hence 
G > {x,y) = G < (y,x). Therefore, p{x,y) = —p(y,x) and 
p{qo,uj q ,T) = —p{—qQ,uj q ,T). The normalisation of po is 
the same as in equilibrium, i.e, for the parametrisation 
(ft) we have 



1 f+°° 

7T~ / qoPo(qo,u q ,r) 



dpo(Vq,t,t') 



dt 



-1 (20) 



It is not difficult to check (^), by using for instance 
© and Po 9 K,f,t) = p (u> q> t~t) = 0. 



III. NEXT TO LEADING ORDER PROPAGATOR 

We will now obtain the NLO correction to the two- 
point Green function (p"3|). For that purpose, we need 
the action in (0) up to four-pion terms: 



S[tt] = S [ir] 



dh 



c 



1 



f 2 (t) 



h2,2 (fit) f 2 (t) 

2 { ' \f(t) P{t) 



o(n e ) 



(21) 



The NLO correction to the two-point function is given 
by the tadpole diagram a) in Fig.|^. It is important to 
bear in mind that in principle we should also include 
a tree diagram coming from the lagrangian, which 
would take care of renormalisation. We shall ignore that 
contribution here and we will justify this approximation 
in section IV in the context of the short-time regime. 



b) 



FIG. 2. One-loop diagrams contributing to a) The NLO 
pion propagator and b) The NLO axial-axial correlator 

After summing over all isospin indices in the loop, 
we can use the differential equation (pT|). However, in 
doing so we have to deal with the divergent quantity 
<5( 4 )(0). We will use dimensional regularisation, which 
is a suitable scheme to deal with ChPT j|, so that 
(5( d )(0) = 0. Therefore, we can replace O z Gq(z,z) — 
—m 2 (z°)Go(z,z). Let us concentrate on t,t' £ C% in 
Fig.|l| (i.e, the Gn component of the propagator [^U) and 
take t and t' positive, which is actually the relevant case 
for our purposes, since the nonequilibrium effects start 
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for positive times. Then, in the mixed representation we 
obtain 

^JVXO,ll(*>* ) = ^0,11 ) 



2u n 



■ coth 



G eq 

-^cos [oJ q (t + t')} 



/ 2 (0) 



+ i 



d?[Ai(t,w,)G^(«,t)G^(t;i / ) 



+ A 2 ®G>(t,i)G>(i,f') 

[A x (t, o; g )G^ (f , QG^ (t, t') 
f A 2 (i)G<(i,i)G<(M') 

di[A 1 (t, Wg )G < (f,t)G >(t,i') 



+ A 2 (t)G<(t,t)G>(t,t')]} 



(22) 



and G^ LO n (i, t') = G^ LO n(t', t), where we have sup- 
pressed for simplicity the uj q dependence of the propaga- 
tors, the dot denotes d/dt, 



Fit) 



'lit) 



m 
.nt) 



at). 



uj z q ) G (t)-2G (t)+4^G (t) 



fit) 



Mt) 



Go® 
Pit) 



(23) 
(24) 



Gq(z°) = Gq(z,z) is the time-dependent equal-time cor- 
relation function and Gg 9 = Go(£ < 0), which is t inde- 
pendent and is given in dimensional rcgularisation by 



u 



jd-l 



q 1 



(2?r) d - 1 2lu c 



coth 
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(25) 



The equilibrium result is finite in the chiral limit, 
where there is no need to rcnormalisc the tadpole di- 
agrams in dimensional regularisation |^]. However, 
out of equilibrium, Go(t) is in general a divergent time- 
dependent quantity, even in the chiral limit, which means 
that the infinities have to be absorbed in the new 
noncquilibrium counterterms we discussed in section p. 

To derive (^2|) we have used the equilibrium result Qlq ) 
to evaluate explicitly the t < integrals. We observe 
that ( p2] ) is ti and tf independent, which is a good check 
of consistency of our result. It can be checked that the 
result is also ti and tf independent for G12, G21 and G 22 . 
Another consistency check is obtained by replacing the 
equilibrium Go(ui q , t, t') and Go(t) — ► Gq 9 in (E2), so that 
we recover the equilibrium result 



^NLO.iiit - — Gl q n(t - t') ^1 



T 2 
12/ 2 



(26) 



The above result agrees with S. Notice that, in turn, 
we have derived the correct equilibrium result using the 
general contour G, which includes both imaginary-time 
and real-time thermal field theory. 



IV. THE NONEQUILIBRIUM PION DECAY 
FUNCTIONS 

At T = the pion decay constant is customarily de- 
fined through the PCAC low-energy theorem for the ex- 
pectation value of the axial current between the vacuum 
and an asymptotic pion state and one relates that am- 
plitude with Green functions via the LSZ formula [|. 
However, in a thermal bath, the concepts of LSZ and 
asymptotic states are subtle and it is much more con- 
venient to work directly with thermal Green functions. 
The object of interest for our purposes is the axial-axial 
correlator 



A<»{x,v) =<T c A-(x)A b ti (y)> 



(27) 



At T = 0, the longitudinal part of the axial-axial spec- 
tral function is proportional, in the chiral limit, to a 
(5-function, signaling the existence of massless NGB (see 
below) and the proportionality constant is f n . However, 
at T 7^ that definition is more subtle: firstly, because of 
the loss of covariance, the tensorial structure compatible 
with current conservation is more general than at T = 
ptjf and one can define two independent and complex 
(spatial) and /* (temporal) Secondly, the pion 
dispersion relation is not in general a (5-function as in 
vacuum. Nevertheless, a proper definition of fn(T) can 
still be given S, which, to one- loop gives the result (|l|) 
with /« = /* = / w real. To higher orders, /* ^ /* and 
Im/^ , Im/£ are nonzero 0] . 

Let us then analyze the axial-axial correlator ( p7j ) in 
our non-equilibrium model. As we said above, the rel- 
evant quantity, as far as / w is concerned, is not (j2^|), 
but its spectral function p a ^ v {x,y), defined by subtract- 
ing its advanced and retarded parts. Again, wc write 
Put = "Pttv by isospin invariance and we concentrate 
on the fast Fourier transform p (UI /(g , o, q, t). Firstly, we 
will write down the most general form for compatible 
with the symmetries and analyse the restrictions imposed 
by the axial current conservation Ward Identity (WI). 

From the definition of the spectral function, we readily 
realise that p^ v (qo,q,T) = -^(-go, -<?> ">")• Then, from 
rotational invariance, 

p tJ (q ,q,T) = qiqjPL(qo,u) q ,r) + 5ijPd(qo,Uq,T) (28) 

with pl.Aqo) = —pL,d(—qo) B On the other hand, p jo is 
a 3-vector, so that 



1 pL and pd correspond in equilibrium and in the notation of 
i, to sgn{q°)p\ql/u>lq 2 and sgn(g°)p3i respectively. 



G 



(29) 



and pjo(qo,q,T) = qjp s (~qo,u q ,T). Thus, p^ v is char- 
acterized, in principle by the four functions pl, Pd, Ps 
and poo- However, not all of them arc independent. 
We have seen that the axial current is classically con- 
served in our model and, therefore, the WI d^p t _ iu (x, y) = 
dvpfivixji/) = must hold (there are no axial flavour 
anomalies for Nf = 2). It can be checked, following 
standard path integral methods Q, that the same equi- 
librium WI holds out of equilibrium. Using the symmetry 
properties discussed before, the WI yields 



q°p 00 - uj^ps - -poo = 
q°ps - uJqPl + t;Ps + Pd = 



(30) 



where the dot denotes d/dr. Thus, only two components 
of Pyn, are independent, as in equilibrium |2^,||. Notice 
that the time derivatives are not present in equilibrium. 

At T = one has p L = 2irf%sgn(q°)5(q 2 ), which de- 
fines /tt and states that in the chiral limit there are NGB 
with the same quantum numbers as the axial current, 
according to Goldstone theorem. That is not the case 
in equilibrium at T ^ p6|,^|, where to define properly 
f-n{T) requires taking the uj q — > + limit, in which Gold- 
stone Theorem requires that a zero energy excitation still 
exists H . Extending the above ideas to non-equilibrium, 
we will define the spatial PDF as 



2 1 f°° 

= 7T lim / dq q pL(qo,u,t) 
2vr^o+ J_ 00 



lim i —p L (u;,t,t') 
ui^a+ at 



(31) 



where in the second line we have written pl in the mixed 
representation. In equilibrium, the above gives / w inde- 
pendent of lu and then there is no need of taking the 
lu — > + limit pl, meaning that the pion distribution 
function still behaves as a <5-function to this order of ap- 
proximation. It is not obvious in the least that the same 
has to hold out of equilibrium to one-loop, as it is indeed 
the case (see below). 

Let us check the consistency of the definition ([Hj). To 
leading order, from (ph, we have 



pi u (uj g ,t,t')=if(t)f(t')p (u; q ,t,t') 



(32) 



with po the spectral density ( |l9| ) to leading order. Then, 
using (|o|), we readily find f°{t)\ = f 2 (t), i.e, the PDF 
coincides with f(t) to leading order, as it should be. Fol- 
lowing the same ideas, we will introduce 



1 f°° 
#(*)/*(*) = 5- I™ / dq ps(qo,u,t) 

27T J_ x 



(33) 



lim 

2-7T w^0+ 



dqoqoPs{qo,u,t) (34) 



The function f*(t) is the nonequilibrium extension of 
the spatial pion decay constant, whereas g^if) vanishes in 
equilibrium. However, according to our previous discus- 
sion on the WI, the three PDF we have introduced are 
not independent. In fact, from their above definitions 
and integrating in go m (|30|) , we find 



(35) 



and then there are only two independent PDF, also as in 
equilibrium 

As we did above with f^(t), we can check now the 
above expressions to leading order. From (^J) we find 

fl(t) L o = f°(t) L o = f(t) and g*{t) LO = f(t), so that 
our definitions are consistent to LO. 

To NLO, we need the axial current to 0(ir 3 ). From (||) 
we have 

A^(x,t) - -f{t)d^ a + <W(*K " ^ (* a d^ 2 

] f(t) 



n 2 d^n a - ^qMttV ] + 0(n 5 ) (36) 



Therefore, according to the chiral power counting ex- 
plained in previous sections, we have two types of con- 
tribution to the axial-axial correlator NLO corrections. 
The first one comes from the NLO correction to the pion 
propagator, which we have evaluated in section [II, in the 
product of the 0(tt) terms above, whereas the second is 
the product of the 0(ir) with the (D(ir 3 ). These two con- 
tributions are represented in Fig.|^ by the diagrams a) and 
b) respectively. As we did in the evaluation of the NLO 
propagator, we are not considering the contributions 
(see below). After evaluating the loops, we obtain p^ v to 
NLO. We give here the result for the spatial-spatial com- 
ponent: 



NLO I 



K, t, t') = p L L {uj qi t, t') + lf(t)f(t')p tad {cj q , t, t') 



fit) r (fl] 



fit') 



fit) 



G (t) 



(37) 



with pl° in (|32j) and p tad = G^ LO - G^ LO - p . The 
Ptad contribution in the above equation comes from the 
diagram a) in FigHand the rest from diagram b). From 
our definitions (plj) and ( |33| ) we obtain 



2 

NLO 



fit) 



i[G (t)~ fimit)} 



im)gAt)} 



NLO 



fit)f\t) 
Jit) 



(38) 



(39) 



with 
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i>— >0+ OX 



G (t) 
Pit) 



(40) 



where we have made use of (p2|), taking, without loss 
of generality, both t, t' S Ci and positive (we know the 
equilibrium answer for the negative real axis). We realise 
that Ai (t, ui) in ( p2|) does not contribute to H(t) and that 
the whole NLO result for the PDF is independent of u>, so 
that there is no need of taking the lo — > + limit, which 
together with f%(t) = flit) in j3q ) and according to our 
previous discussion, we may interpret as the existence 
(to this order) of undamped NGB in the nonequilibrium 
plasma, propagating at the speed of light (see Q ) . 

In addition, when we substitute ( ft0| ) in ( p8| ) and (39), 
we see that g v (t) to NLO satisfies the WI relation (35), 
which is a consistency check. We finally obtain 

\fW)]%LO = [fMlLO = - 2 * G °W ( 41 ) 

This is the one-loop relationship between the PDF and 
fit) extending the equilibrium result (|l|), which we 
recover (for the general contour C) simply by replac- 
ing (|25|) in Q). We want to stress that the PDF are 
observable and therefore (^lj) should be independent of 
the parametrisation chosen for the pion fields. We have 
checked it explicitly by calculating with the tt fields de- 
fined before. Both the propagator and the axial-axial 
spectral functions change, but (|4l]) remains the same 
with the same definitions (|3|) and ([33]) and the same 
Go(t). In turn, let us mention that by extending naively 
the low T results in the LSM (see ||) to nonequilibrium 
we would have v 2 (t) = f%(i) - iG (i) with v(t) = (a)(t). 
We will make use of this assumption below to compare 
our results with p6| . 

Notice also that it is licit to replace f(t) by f n (t) 
for Go(t) in (f4l|), to the order we are considering here. 
Therefore, using (|4l]) allows to express all physical non- 
equilibrium observables (like decay rates, masses, etc) 
obtainable from our model to one loop in ChPT, in 
terms of the physical f w (t), which could be determined, 
for instance, by analysing non-equilibrium lepton decays 
7r — > Ivi. A very important aspect in this program is 
renormalisation. As we said before, Go(t) contains in 
general time-dependent infinities. However, as we are go- 
ing to see, in certain regimes (for instance, at short times) 
it is still possible to find a finite answer for Go(t) in di- 
mensional regularisation (as in equilibrium in the chiral 
limit) so that the counterterms coming from are not 
needed. There still could be finite contributions coming 
from £( 4 ) but those are typically of 0(1O~ 3 ) with respect 
to those coming from the loops so that our ap- 

proximation of neglecting would be justified in that 
case, which is the one we will consider in the remaining 



of this work. Explicit knowledge of f(t) would require 
to solve the hydrodynamic equations for the plasma, or, 
equivalently, the Einstein equations for the metric. Alter- 
natively, one can follow the approach of treating f(t) as 
external, choosingfor instance different parametrisations 
compatible with (g). This is left for future investigation. 

For short times, the particular form of f(t) is not im- 
portant and we can parametrise the nonequilibrium dy- 
namics in terms of a few parameters, which are the val- 
ues of fit) and its derivatives at t = 0. This is equivalent 
to describe the expansion of the Universe in terms of the 
Hubble parameter, the deceleration parameter and so on. 
As we discussed in section [n], this approach is justified 
for times t < t max with t max ~ l//„.(0) ~ 2 fm/c, which 
are indeed not so short, regarding the typical time scales 
involved in a RHIC (see our comments below and in sec- 
tion |) . Let us then start by considering the differential 
equation ( p"7j ) with KMS conditions at tj. In addition, 
as we have discussed in section [nj we shall impose that 
Go(uj q ,t,t') be continuous and differentiable in t — and 
t' = 0. The solution of (17) can be written in terms of 
two arbitrary solutions <7j(i) of the homogeneous equa- 
tion, satisfying the Wronskian condition <?i<72 — — L 
In pfj, the general solution of ( |l7| ) satisfying KMS con- 
ditions was constructed in terms of the g%{$) functions, 
for t and t' in the contour C (as we have said before, 
it is enough for us to take t, t' <E C\). We remark that 
the solution of the homogeneous equation is only known 
explicitly for a some special choices of m 2 (t) |^0|,^2) . 

Therefore, we expand both fit) and the solution of 
the homogeneous equation near t — 0. Both the dif- 
ferential equation and the Wronskian condition reduce 
then to algebraic relations between the first coefficients 
in the expansion. With the boundary, continuity and dif- 
ferentiability conditions discussed above, the solution is 
uniquely defined and we find to lowest order 



G (u)q,t,t') = -- — coth 

ZLOq 



[l-m 2 t 2 +0(mH 4 )] 
(42) 



2 Notice that Go it) in ( ftl] ) depends implicitly on fit), 
through Q. 



with m 2 = m 2 (0) = -/(0)//(0). For m 2 < we see 
the onset of unstable modes, making the pion correla- 
tion function grow with time. We also observe that for 
short times, the time dependence of the mixed propaga- 
tor factorises. Therefore, the momentum dependence is 
the same as in equilibrium and then we can integrate it 
in dimensional regularisation, using (5q), to get a finite 
answer. As we said before, in this regime we are able to 
avoid renormalisation. Then, from (|4l]) and ( |42| ) we get 

flit) = f 2 [ain) + 2Ht - c(T,)i 2 + Oi\m\H 3 )] (43) 

where / = /(0), H = /(0)//(0) (H and m are C?(p) and 
play the role of the Hubble and deceleration parameters 
in the expansion, respectively) and 

T 2 

aft) = 1 - ^ 
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c{T % ) = m 2 a(^) - H 2 



(44) 



where T c = V&U - 228 MeV. Therefore, @ pro- 
vides the noncquilibrium relationship between the LO 
Taylor coefficients, independent of T, appearing in the 
lagrangian and the physical IVdependcnt NLO ones, ex- 
tending the equilibrium formula when / is the only 
lagrangian parameter. The coefficients of order zero and 
two in the Taylor expansion of / w (t) are corrected by the 
same one-loop X^-dependent factor, whereas the first or- 
der coefficient H does not acquire NLO corrections. We 
would like to remark that the corrections coming from 
£W that we have neglected are, in addition, independent 
of Ti, since they only contribute at tree level to NLO. 

We will proceed now to discuss some physical effects 
related to the behaviour of f n (t). For that purpose, and 
based upon |l]), let us define a time-dependent effective 
temperature T(t) as 



i 



T 2 {t) 



(45) 



i.e, we parametrise the nonequilibrium effects of fn{t) in 
T(t). Let us define in addition the critical time t c , as 
MQ = (T{t c ) = T c ) and the freezing time T{t f ) = 0, 
fir iff) = f ■ These are the relevant time scales when the 
system is heating or cooling down, respectively. The sign 
of H determines whether the plasma is initially heated 
(H < 0) or cooled down [H > 0). However, it is the sign 
of c(Tj) in (|4i|) what determines the behaviour at longer 
times. Thus, we will follow the short-time evolution of 
f-n{t) until it reaches either t c or tf, imposing that < 
T(t) < T c . We are following a similar approach as that 
of equilibrium when one extrapolates ([!]) until T c . 

In order to estimate the above time scales, let us take 
H 2 ~ \m 2 \ and retain only the leading order in x = 
T 2 /T 2 , consistently with the chiral expansion. Then, we 
have the following cases: 

1. H > 0: Cooling down until it ~ x/2. 

2. H < 0, m 2 > 0: Heating until t c ~ (1 - 3af/4)/2. 

3. H < 0, m 2 < 0: Heating until i m ~ (I + x/2)/2 
and then cooling down until if ~ l + x, passing 
through the equilibrium time t eq = 2t m such that 

MM) = MO)- 

4. H = 0, m 2 > 0: Heating until t 2 m 2 = 1, indepen- 
dent of Ti. 

5. H = 0, m 2 < 0: Cooling until t 2 m 2 ~ — x(l + x) 

where i — t\H\. We observe that the effect of the unsta- 
ble modes (m 2 < 0) is always to cool down the system. 
The freezing time tt for H < is much longer (accord- 
ing to our power counting) than that for H > 0, which 
is 0(x/H). Indeed, the 0(1/H) time scales are much 
longer than those to which our short-time approxima- 
tion remains valid and they have to be understood as 



estimates, similarly to the equilibrium case, where one 
estimates the critical temperature, even though the low 
T approach breaks down for T ~ T c . Taking typical val- 
ues T, ~ \H\ ~ \m\ ~ 100 MeV (x ~ 0.16) we get tf ~ 
0.2 fm/c for cases 1 and 5, t c ~ 2 fm/c for cases 2 and 4 
and i/ — 2.3 fm/c for case 3. 

On the other hand, evaluating the above time scales 
for initial temperatures closer to T c , by expanding now 
for small a(Ti ) in @), we get if ~ V2 - 1 for case 1, 
f c ~ a{Ti)/2 for case 2, t OT ~ 1, ~ l + \/2 for case 3 and 
t 2 m 2 = —x/a(Ti) for case 5. We can compare with v(t) in 
p6[ (see our comment above), using their initial values 
Ti ~ 200 MeV and ~ 400 MeV (which are clearly 
too high for our low-energy approach). We see that our 
estimates for the time evolution duration are of the same 
order as those in p6| , although ours are somewhat lower. 
An important remark is that in typical LSM simulations 
like [lql , v(t) reaches a stationary value, about which it 
oscillates, the relevant time scale being the time it takes 
the system to reach such value. It is clear that we cannot 
predict that type of behaviour only within our short-time 
approach, quadratic in time, but only estimate the time 
scales involved [J i.e, we cannot see the system thermalise 
in that sense. 



V. CONCLUSIONS AND OUTLOOK 

We have extended chiral lagrangians and ChPT out of 
thermal equilibrium. The chiral power counting requires 
all time derivatives to be 0(p) and to lowest order our 
model is a NLSM with / — * fit). This model accommo- 
dates unstable pion modes and corresponds to a spatially 
flat RW metric in conformal time with a(t) = /(£)// (0) 
and minimal coupling. For a non-minimal coupling, the 
chiral symmetry would be broken and conformal invari- 
ance restored for £ = 1/6. 

We have applied our model to study the time depen- 
dent PDF, which extend the equilibrium pion decay con- 
stants. We find that in general there are two indepen- 
dent PDF, as in equilibrium. To one loop in ChPT such 
functions coincide and we have given their expression in 
terms of the equal-time correlation function, analysing 
the lowest order short-time coefficients of We have 

studied the relevant time scales and their dependence 
with T, paying special attention to the role of the unsta- 
ble modes in cooling down the plasma. In this work we 
have treated nonequilibrium renormalisation by remain- 
ing within the short-time expansion, where the equal- 
time correlator turns out to be free from UV divergences 
in the chiral limit in dimensional regularisation and hence 
it is reasonable to neglect the contributions coming from 



3 Similarly as to why ChPT cannot see the phase transition 
(see comments in section tth 
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£( 4 ) in our analysis. Nonetheless, we have seen that, 
thanks to the analogy with curved space-time, we have a 
well-dehned procedure to construct so as to absorb 
all the loop divergences, which in general will be time- 
dependent, and find a finite answer for the observables. 

Among the aspects we would like to study further is 
to include the nonequilibrium tree contributions, so 
that we can renormalise fn(i) Vt and analyse the long 
time behaviour properly. An interesting remark is that 
for long times with unstable modes of negative mass 
squared m 2 , only the k 2 < \m 2 \ modes are important, 
so that there is an effective ultraviolet cutoff and we can 
forget about UV divergences in that case as well. An- 
other point we leave for future investigation is the be- 
haviour of the two-point correlation function at different 
space points for different choices of f(t), which would al- 
low us to investigate the formation of regions of unstable 
vacua (DCC) in our model. Other applications and ex- 
tensions of our approach, to be explored in the future, 
include photon production in the pion sector (by gaug- 
ing the theory and including 7r° decay photon produc- 
tion through the anomalous Wess-Zumino-Witten term), 
the quark condensate time dependence (by considering 
the mass explicit symmetry-breaking terms we have ne- 
glected here) and the Nf = 3 case. Finally, we hope 
that by implementing resummation methods like large 
N in our model, we may be able to study the long-time 
thermalisation process. 
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